EMC effect and jet energy loss in relativistic deuteron-nucleus collisions 
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We investigate tlie influence of modified nuclear parton distribution functions (PDFs) on high- 
Pt hadron production at RHIC and LHC energies using a pQCD-improved parton model. For 
application at RHIC, we focus on the possible contribution of the EMC modification of the nuclear 
PDFs in the x > 0.3 region to the observed suppression of tt" production at pr > 10 GeV/c 
in dAu collisions. We study three different parameterizations of the nuclear PDF modifications 
and find that they give consistent results for RdAuiPr) for neutral pions in the region 10 GeV/c 
< Pt J; 20 GeV/c. We find that the EMC suppression of the parton distributions in the Au nucleus 
does not strongly influence the RdAu for TfO in the pt region where the suppression is observed. Using 
the HKN parameterization, we evaluate systematic errors in the theoretical RdAuiPr) resulting from 
uncertainties in the nuclear PDFs. The measured nuclear modiflcation factor is inconsistent with 
the pQCD model result for pr > 10 GeV/c even when the systematic uncertainties in the nuclear 
PDFs are accounted for. The inclusion of a small final-state energy loss can reduce the discrepancy 
with the data, but we cannot perfectly reproduce the pr dependence of the measured RdAuipr). 
For the LHC, we find that shadowing of the nuclear PDFs produces a large suppression in the yield 
of hadrons with pT < 100 GeV/c in p{d)A collisions. 

PACS numbers: 24.85.+p,25.75.-q, 12.38.Mh 



I. INTRODUCTION 

Deuteron-gold (dAu) collisions have been extensively 
studied at the Relativistic Heavy Ion CoUider (RHIC), 
both for their intrinsic interest and as a control experi- 
ment [l|, [3, 3 to judge the suppression seen in central 
gold-gold (AuAu) collisions at sufficiently high trans- 
verse momenta (pt) Unexpectedly, not only AuAu 
data, but recent extended px coverage dAu data also dis- 
play a suppressed nuclear modification factor in central 
collisions [y]- This motivates a study of possible mecha- 
nisms that may result in a nuclear modification factor 
smaller than unity at sufficiently high transverse mo- 
menta {pt > 6 GcV/c) in central dAu collisions. 

The nuclear modification factor RdAu compares the 
spectra of produced particles in dAu collisions to a hy- 
pothetical scenario in which the nuclear collisions are as- 
sumed to be a superposition of the appropriate number of 
nucleon-nucleon collisions. In the transverse momentum 
window 2 GeV/c ^ Pt ^ 6 GeV/c, the nuclear modifi- 
cation factors are dominated by the Cronin peak 0, Q . 
Several physical pictures of this enhancement have been 



proposed iJM El [ll [H, [3. One family of mod- 
els [HI, [H, 17, 13, S 2C| advances an explanation of the 
Cronin effect in terms of the interplay between nuclear 
shadowing [HI, [H, [H, [H, [H, and the multiple scat- 
tering of particles propagat ing in the strongly-interacting 
medium (multiscattering) [3 [itI- jisl- [ist ■ Using the HI- 
JING shadowing prescription |23l. |24|. our model gave 
a reasonable description of the Cronin effect in central 
collisions at midrapidity [23, . At the same time, we 
obtained nuclear modification factors close to unity at 
high transverse momenta (6 GeV/c < Pt ^ 20 GeV/c). 

In the high-pT region multiscattering no longer affects 
RdAu- It is then natural to ask if nuclear shadowing 
can explain suppression effects at high px- In particular, 
since at ^Jsnn = 200 ^GeV we are in the EMC region of 
the shadowing function for 6 GeV/c Pt ^ 20 GeV/c, 
we ask if the EMC effect [2l| plays a role in understanding 
these experiments. 

In this paper we first investigate dAu collisions at the 
highest RHIC energies and the role of the EMC effect at 
transverse momenta up to 50—70 GeV/ c. We then exam- 
ine whether recent shadowing parameterizations incorpo- 
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rating theoretical uncertainties for the first time[25l. |26| 
can account for the experimental information. We also 
display calculational results with a modest energy loss us- 
ing energy-loss parameters applied earlier to AuAu data. 
Finally we extend our considerations to the energy range 
of the Large Hadron Collider (LHC). 



II. THE EMC EFFECT IN DEUTERON-GOLD 
COLLISIONS 

Figure [U displays recent PHENIX data (triangles with 
error bars)[l] for the most central dAu collisions, where 
a high-pT suppression is clearly seen. While incoher- 
ent multiscattering can only lead to enhancement in 
RdAm nuclear shadowing displays two regions where an 
RdAu < 1 can be expected: (i) at small x {x_< 0.2), 
and (ii) in the EMC region (0.5 < a; < 0.9) j^. At 
RHIC energies the small-x region is inconsequential at 
> 6 GeV/c. Thus we focus attention on the EMC 
effect as a possible mechanism for the measured suppres- 
sion. Various shadowing parameterizations developed in 
the last 15 years [12, ^ [H, ^ show different be- 
haviors at small-x, but the EMC region appears rather 
robust in most models. 

To see the effect of the EMC region, we calculate pion 
production in a wide momentum range. For this pur pos e 
we use a pcrturbative QCD improved parton model [l8| . 
The model is based on the factorization theorem and 
generates the invariant cross section as a convolution of 
(nuclear) parton distribution functions fa/Aj pcrturba- 
tive QCD cross sections dCT''''^'^''/d^, and fragmentation 
functions D^/^. We perform the calculation in leading 



order, following Refs. [18|, M, M, M, M, \M 
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where and represent the factorization and frag- 
mentation scales, respectively, Xa-, Xb, and Zc are mo- 
mentum fractions, and k^-s stand for two-dimensional 
transverse momentum vectors. The initial state effects 
of shadowing and multiscattering are included following 
the treatment in Refs. [isl. [27l. [23] . 

Since the effects we investigate are on the 10 — 20% 
level, it is customary to present the obtained results on 
a linear scale in terms of the nuclear modification factor 
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Here {Nbin) is the average number of binary collisions in 
the various impact-parameter bins. 

Together with the data in central dAu collisions, we 
display our results with several shadowing parameteriza- 
tions in Fig[TJ We use the HIJING shadowing including 
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FIG. 1; (Color online) The nuclear modification factor, RdAu 
in central (0 — 20 %) dAu collisions for tv" . Data are from 
Ref . . Theoretical results are calculated with different shad- 
owing parameterizations (see text). 
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FIG. 2: (Color online) The influence of the uncertainty in the 
HKN shadowing parameterization [2^ on the factor RdAu in 
the most central dAu collisions. Data are from Ref. 0|. 



nuclear multiscattering (solid lines), the EKS shadow- 
ing (where multiscattering is represented by strong anti- 
sliadowing) (dashed line), and the HKN parameteriza- 
tion (with and without nuclear multiscattering, dotted 
and dash-dotted lines, respectively). It can be seen in 
Fig. [T] that the suppression associated with the EMC ef- 
fect shows up at transverse momenta px ^ 20 GeV/c in 
all models considered, and does not explain the suppres- 
sion in the data at around 10 GeV/c. The Cronin peak 
at « 3 GeV/c is best reproduced by the HIJING 
parameterization. The "HKN+multiscattering" model 
appears to overshoot the data at low px, while it gives 
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FIG. 3; (Color online) RdAu for (0-20 %) dAu collisions with 
jet energy loss at opacity L/\g = 1.0. Data are from Ref. 0|. 



very similar results to HIJING at 5 GeV/c. Since 

the uncertainty of the HKN nuclear parton distribution 
functions is also available, we next examine the theoret- 
ical uncertainty of the description. 

In Fig.[2]the effect of the uncertainty given by the HKN 
shadowing parameterization is illustrated on the nuclear 
modification factor, RdAu- The errors are calculated by 
the Hessian method, using the original code of the HKN 
group. The mean value as a function of pt is represented 
by a solid line, surrounded by an error band of approx- 
imately ±10%. Although the data fall within this band 
at low Pt, for pt^^ GeV/c the observed suppression is 
stronger than allowed by this calculation. 

Since we are not aware of any other initial-state sup- 
pression, we consider physics operating in the final state. 
In particular, jet energy loss suggests itself [H, H^. This 
idea is supported by the presence of a minor energy 
loss effect in peripheral AuAu collisions [s^- Assum- 
ing L = 1.5 fm for the average static transverse size of 
the traversed medium, and the usual GLV parameters of 
A = 1.5 fm mean free path and /x = 0.5 GeV screening 
mass, we calculate the effect of jet energy loss on the nu- 
clear modification factor for central dAu collisions. This 
is displayed in Fig. [3] It can be seen that inclusion of 
this energy loss results in a parallel down-shift of RdAu 
relative to curves in Fig O The slope of the data as a 
function of pr is still very different from that of the calcu- 
lated results. Nevertheless, taking into account all exper- 
imental and theoretical uncertainties, one could consider 
the displayed result with L/A = 1 to be an acceptable 
compromise. 

Jet energy loss depends on the transverse parton den- 
sity, which relates to the measurable hadronic quantity 
1/A± ■ dNch/dy, where A± is the transverse area of the 
deconfined region. Assuming a realistic geometry for cen- 
tral AuAu and dAu collisions and considering the exper- 
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FIG. 4: (Color online) The nuclear modification factor for 
pion production in — 10% most central dPb collisions at dif- 
ferent energies. At y^sjvjv = 8.8 TeV the thin lines illustrate 
the effect of increasing energy loss. 

imental data on dNch/dVi we obtain only a factor of 2 
difference for the transverse parton densities between the 
AuAu and dAu cases. On this basis, one could expect 
the jet-quenching effect in dAu to be even stronger than 
shown by the calculated band in Fig. 3. However, jet 
energy loss in non-thermal matter is an open question, 
which we plan to investigate in a forthcoming paper [ssl ] . 



III. PREDICTIONS FOR THE LHC 

We expect that the EMC region will shift to higher and 
higher transverse momenta with increasing collision en- 
ergy. At the same time, due to the increasing parton den- 
sity, the energy loss should also become larger. We repeat 
our calculation for dPb collisions at ^/snn = 900 GeV 
and 8.8 TeV to display these tendencies. The results are 
shown in Fig. 4. Following the ordering of the lines on 
the left of the Figure, the top curve represents the re- 
sults with HIJING shadowing from Fig. [TJ The second 
line corresponds to a similar calculation (i.e. HIJING 
shadowing, no energy loss) at 900 GeV cm. energy. 
The next line shows the result at 8.8 TeV without jet 
quenching, while the dotted and dashed lines illustrate 
the effect of energy loss. First we took L = 1.5 fm for 
the transverse size of the medium, and used the value of 
A = 1.5 fm as earlier. The introduction of a smaller A 
value to represent the increase in the transverse density 
of colored scattering centers with increasing cm. energy 
gives L/X = 3.0. It can be seen that the dip correspond- 
ing to the EMC effect shifts to higher transverse momenta 
with increasing energy as expected. While RdAu is above 
unity at ~ 3 GeV/c for 0.2 TeV (Cronin effect), at 
0.9 and 8.8 TeV, where we are deep in the shadowing re- 
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gion at these transverse momenta, at most a small ripple 
appears on RdAm which is rising towards one with in- 
creasing transverse momentum in this region. The effect 
of varying A at 8.8 TeV is minimal at transverse mo- 
menta pt ^ 100 GcV/c, but quite significant at 5 GeV/c 
< PT < 10 GeV/c. 

IV. CONCLUSION 

In conclusion, while this solution appears initially 
tempting, the EMC effect docs not explain the unex- 
pected suppression seen at ~ 10 GeV/c in RdAu at 
^/snn = 200 AGeV. Taking into account experimental 
uncertainties as well as the uncertainties of the HKN nu- 
clear parton distribution functions, the experimental and 
theoretical results can be brought into agreement using 
a non-negligible amount of final state parton energy loss 
(with standard opacity parameters). The non-thermal 
nature of the dAu system casts some doubt on the pa- 
rameter values. The presence of jet quenching in the dAu 
system is somewhat surprising at first sight, but is jus- 
tified by the similarity (within a factor of 2) of the real 



transverse densities in the AuAu and dAu systems. 

For LHC energies we do not have a fully reliable base- 
line calculation at present. It is clearly seen, however, 
that shadowing (i.e suppression) will dominate the mo- 
mentum region up to px S 100 — 200 GcV/c. The effect 
of final state interactions (jet energy loss) may also re- 
duce the nuclear modification factor, especially in this 
momentum region. The EMC dip in the nuclear modifi- 
cation factor moves toward higher and higher transverse 
momenta with increasing energy. 
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